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I.     Introduction. 

In  a  previous  paper  [1],  we  reported  on  a  simple 
experimental  algebraic  manipulator  which  consisted  of  a 
set  of  subroutines  for  doing  elementary  algebra  with 
expressions  which  were  essentially  polynomials.   The 
information  gained  from  experience  with  that  system,  and 
from  reports  of  other  workers  in  the  field  [3,4],  led  us 
to  design  an  enlarged  system  which  would  include  arbitrary 
fixnctions  of  one  variable  and  differentiation  capabilities 
as  well  as  a  more  problem-oriented  appearance. 

The  original  program  (Symbolang)  was  rewritten 
using  a  revised  internal  representation  of  expressions. 
The  new  package  forms  the  basis  for  the  current  system, 
called  SYMB66.   These  routines  still  use  a  subset  of 
the  list  processor  SLIP,  operate  primarily  on  Fortran-like 
expressions,  and  allow  various  forms  of  symbolic  algebraic 
manipulation.   Superimposed  on  these  routines  is  a  pre- 
processor, so  that  a  "non-expert"  problem  oriented  user 
might  not  have  to  be  concerned  with  the  internal  machina- 
tions of  the  list  processor.   This  allows  the  use  of  a 
more  natural  problem  description. 

Those  interested  in  using  this  higher  level  problem 
oriented  language  will  find  a  telegraphic  description  of  it 
in  Section  II,  and  more  complete  details  in  Section  V. 
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Sections  III  and  IV  deal  with  specific  details  of  the 
system.   Section  III  discusses  machine  level  notation, 
input  and  output,  and  Section  IV  discusses  the  algebraic 
computation  package  of  routines . 
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II.   Instant  SYMB66 


A   preprocessor  is  available  for  SYMB66  which 
relieves  the  user  of  the  task  of  learning  the  internal 
structure  and  other  minutae  of  the  system.  A   summary 
of  simplified  commands  is  given  followed  by  a  detailed 
explanation  of  the  usage  and  function  of  the  processes 
Section  V  gives  more  comprehensive  information  about 
the  preprocessor. 

Comment 


COMMAND 

Abbrev. 
(if  any) 

RESERVE 



EXPRKSSI0N 

EXP 

NAME 



THE  DERIVATIVE  F^ 

TDF 

ADD 



SUBTRACT 

SUB 

DIVIDE 

DIV 

MULTIPLY 

MUL 

DIFFERENTIATE 

DIF 

EVALUATE 

EVA 

GET    C EFFICIENT... 

Gcgl 

GET   PPWER.  .  . 

GP^ 

TRUNCATE 

TRU 

S^VE 



SUBSTITUTE 

SUB 

REPLACE 

REP 

^2IUTPUT 

gfUT 

EMPTY 



DESTR0Y 

DES 

Ci^'PY 



Reserves  space 
Creates  Expression 
Creates  null  expressions 
Defines  derivatives 
Adds  expressions 
Subtracts  expressions 
Divides  by  a  variable 
Multiplies  expressions 
Differentiates  expressions 
Substitutes  constants 
of  variables  to  powers 
of  specified  variable 
Truncates  series 
Solves  simple  linear  eq. 
See  detailed  description, 
another  type  of  substitution 
Prints  expressions 
Clears  expression 
Annihilates  expressions 
Makes  copies  of  expression 
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SIMPLIFY 

IF 

AND  PLACE  (2(lT 

AND  TRUNCATE  0 

EXCLUDING  FUNGTIQlT 

ARGUMENTS 

WITH   RESPECT   Tgi 


SIM 

AFgi 
ATgi 


EFA 

WRT 


Collects  terms 

Used  in  tests 

Used  in  other  commands 

Another  substitution 

option 

Used  in  substitution 
Used  in  DIF  command. 


RESERVE  n  -  Allocates  n  memory  locations  for  symbolic 

manipulation.   Must  be  used  first. 

EXPRESSION  ia  =  "Fortran  type  expression"  (use  of  " /" 

prohitibed)  -  Creates  a  symbolic  expression  £a    . 

NAME  ia,ib,ic,...  (up  to  5  var)  -  Designates  names  of 

expressions  which  have  the  value  zero  until  other  terms 

are  added  to  them. 

THE  DERIVATIVE  FOR  "function  name",  IS  ia  ^  "expression 

for  derivative".  -  Defines  the  derivative  of  the  function 

named  to  be  the  expression  for  derivative  which  is  to  be 

an  expression  with  the  convention  that  the  argiiment  of 

the  function  is  replaced  by  zero  in  the  expression. 

ADD  ia,    TO  i,h   -   Adds  the  expression  £a   to  ib  and 

destroys  la. 

SUBTRACT  £a,    FROM  jb  -  Subtracts  ia  from  ib  and 

destroys  ia. 

DIVIDE  iaj  BY  char^  AND  PLACE  0^    ib  -  Divides  an  expression 

by  a  variable. 
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MULTIPLY  Ha,    BY  Vo ,    APQ'  Ho.   -   Multiplies  the  expression 
ia.   by  the  expression  Hh   and  adds  the  results  to  ^c. 
MULTIPLY  iaj  BY  con  -  Multiplies  the  expression  Ha   by  the 
constant  con  and  leaves  the  result  on  Ha.       £a   is  the 
original  and  resulting  expression. 

The  last  five  instructions  may  be  written  in  the 
form 

^  =  ib  +  ^a     (for  ADD,  SUB) 

ib  =  fe  /  char    (for  DIV) 

ic  =    Ha  *    Kb  (for  MUL) 

Ha   =    Ha  *    con     (for  Mult  by  const) 
provided  that  an  S  appears  in  columoi  1. 

DIFFERENTIATE  Ha,    WITH  RESPECT  TO  char,  APjZT  ^b  - 

Differentiates  an  expression  Ha   with  respect  to  the 

variable  char  and  places  the  result  on  the  expression  ^b. 

EVALUATE  Ha,    F^  char  =  X,    char  =  X,  ...  (up  to  3)  - 

Substitutes  the  indicated  numbers  X  for  the  indicated 

variables  char  on  the  expression  Ha   and  leaves  the 

result   on  Ha. 

GET  G^FFICIENT  gfF   char  **    e,    IN  Ha,    APg^  Hh   -   This  is 

used  to  get  the  coefficient  of  a  variable  (char)  to  a  power. 

If  **e   is  omitted  the  exponent  is  assumed  to  be  1. 

GET  P0WER  0F  char,  IN  Ha,    AP0  var  -  This  is  used  to  find 

the  power  of  a  variable  char  in  the  expression  Ha.      The 

result  is  stored  in  "var."  Ordinarily  char  would  appear 
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in  only  one  term  of  the  expression  £a. 

TRUNCATE  isj  Fgl!R  char  **  e  -  Eliminates  terms  in  an 

expression  ia  that  contain  char  to  a  power  greater  than 

e.   ia  is  the  original  and  resulting  expression. 

SOLVE  £a,    FOR  char,  APO  ib  -  Solves  the  equation  ia  =  0 

assuming  that  char  appears  linearly  in  only  one  term  of 

the  expression. 

The  next  six  explanations  pertain  to  substitution. 
FffiPLACE  is  used  to  mean  insert  an  expression  into  the 
original  in  place  of  a  designated  character.  SUBSTITUTE 
is  used  to  mean  insert  the  expression  and  then  carry  out 
the  necessary  multiplications. 

SUBSTITUTE  £a,    FOR  char,  IN  ih,  APO  ic,  EXCLUDING 

FUNCTION  ARGUMENTS  -  Substitutes  the  expression  ia  for 

char  where  it  appears  in  £b   except  where  it  appears  in 

the  argument  of  a  function. 

SUB  ia,  Fg^R  char,  IN  Zb,    AFgi  £c   -   Substitutes  the 

expression  ia  for  char  wherever  it  appears  in  ibj 

even  in  the  argument  of  a  function. 

SUB   ia,    FOR  char,    IN   ib,    APgf  ic,    EFA  AND  TRUNCATE  ON 

var,  pwr,  var,  pwr,  ...  (up  to  5  pairs)  -  Substituees 

the  expression  ia  for  char  where  it  appears  in  ib_, 

excluding  appearances  in  function  arguments,  and  then 

eliminates  all  terms  containing  the  variable  var  to  a 

power  greater  than  that  specified.   The  results  are  put  on  ic, 
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SUB  ia,  FOR  char,  IN  ib,  APg^  ^c,  ATg^  var, power, .. . 

(up  to  5  pairs)  -  Substitutes  the  expression  ia  for 

the  variable  char  wherever  it  appears  in  ib  and  puts 

the  results  on  2,c   after  truncating  terms  in  which  var 

appears  to  a  power  greater  than  that  specified  in  the 

sequence  var,  power,  var,  power,  ...  . 

REPLACE  char,  IN  ia,  WITH  ib,  APg^  ic  -  Simply  replaces 

char  wherever  it  appears  in  ia  by  the  expression  ib. 

The  results  are  put  on  ic. 

REP  char,  IN  ia,  WITH  ib,  APg^  ic,  EFA  -  Replaces  cher 

with  ib  where  it  appears  in  ia  except  where  it  appears 

in  a  function  argument.   The  results  are  added  to  ic. 

OUTPUT  ia , ib , —   -  Prints  one  to  five  expressions  ia, ib, 

...   in  Fortran-like  format. 

EMPTY  ia,ib,...  (up  to  5)  -  Eliminates  all  terms  of  an 

expression  making  the  expression  identically  0. 

DESTROY  ia,ib,...  (up  to  5)  -  Erases  the  expressions 

from  memory. 

COPY  ia,  Tg^  ib  -  Makes  a  copy  ib  of  the  expression  ia. 

SIMPLIFY  ia, ib,...  (up  to  5)  -  Collects  terms  in  an 

expression  including  function  arguments. 

IF  ia,  C^STANT  var  =  CON  "statement  number"  -  If  the 

expression  ia  is  a  constant,  the  value  of  the  constant 

is  stored  in  var  and  control  transferred  to  the  indicated 

statement  number.   Otherwise  control  is  transferred  to 
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the  next  statement. 

IF  £a,    ZERgI   statement  number  -  If  la   is  equal  to  zero, 

control  is  transferred  to  the  statement  indicated. 

Otherwise  control  is  transferred  to  the  next  statement. 

IF  £a,    EQUAL  £b,    statement  niimber  -  If  £a   is  an 

exact  copy  of  ib,  control  goes  to  the  statement  number 

specified;  otherwise  control  goes  to  the  next  statement. 

A  Simple  Example: 

RESERVE  3500 

PRINT  10 
10   FORMAT (*  THIS  IS  THE  BEGINNING  OF  THE  RUN*^  ) 

NAME  SINY,  LTgfZ 

EXP  SINY  =  Y  -.166  *  Y  **3  +  .  0085555*  Y  **  5 

EXPLTi2l^  =  l+Z  +   .5*    Z**2+.l66*    Z  **    5 

SUBSTITUTE  SINY,    Fg(R  Z,    IN   LTgfZ,    APgf  RESULT 

PUTPUT   RESULT 

CALL  EXIT 

END 
The  computer  output  is  I 
THIS  IS  THE  BEGINNING  ^  THE  RUN 

SINY  =  Y  -  .166  *Y  **5  +  .0085533  *Y  **5 

LT^  =  l  +  Z  +  .5*  Z**2  +  .166*Z**5 
RESULT  =1.  +  Y  +  .500000000  *Y  **2  -  .I66OOOOOO  *Y  **4. 
-.074554700  *Y  **5  +  .022111500  *Y  **6.    +  .017872871  *Y  **Y. 
-.001385528  *y  **8  -  .002157128  *Y  **9.  +  .000054722  *Y  **10. 
+.000148940  *Y  **11-  .000005741  *Y  **15.+  . OOOOOOO96  *Y  **15. 
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Notes  on  Example: 

Line         Comment 

1  reserves  space  and  sets  up  COMMON 

2,3  FORTRAN 

4  creates  null  expressiorE  named  SINY  and  L(2i^Z 

(not  necessary) 
5:,6         creates  expressions  --  an  alternative  way 

would  be  to  omit  EXP  and  insert  X  in  col.  1 

7  APO  means  "and  put  on";  RESULT  is  created 
explicitly;  SUBSTITUTE  may  be  abbreviated 
as  SUB 

8  Prints  expression  RESULT 
9,10        FORTRAN 

Since  division  is  not  allowed  in  input,  a   subexpression 
which  occurs  as  a  divisor  should  be  written  with  an 
exponent  of  -1. 
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III.   Notation. 

A.      Concept  of  an  Expression. 

From  the  user's  standpoint,  the  basic  unit  of 
operation  is  the  "expression"  which  is  composed  of  one 
or  more  terms.   While  the  collection  of  algebra  routines 
is  organized  primarily  to  compute  with  these  expressions, 
and  the  programmer  generally  need  not  concern  himself 
with  individual  terms,  the  user  may  have  interest  in 
operating  on  specific  terms  and  thus  an  explanation  of 
the  structure  of  expressions  as  they  appear  on  SLIP 
lists  [2]  follows. 

Representation  of  an  Expression  on  Lists. 

The  algebraic  manipulator  is  designed  to  deal  with 
those  expressions  whose  Fortran  representation  contains 
neither  division  sign  nor  variable  exponentiation  nor 
function  with  more  than  one  argument.   An  example  of  such 
an  expression  is 

(1)   1.5  y?'^   Y  Sin  (X^)  Sin^'5(Y^)  Cos  (X+Y- Sin(Z) )+l . 2X^Z 

These  expressions  are  contained  on  lists  of  the  SLIP 

list  processor  [2]  in  which  data  which  are  pointers  to  lists 

are  immediately  identifiable  as  such.  The  representation 
is  as  follows. 
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Expression  -  A  list  of  pointers  to  sublists. 

An  expression  is  a  siim  of  the  terms  represented  on 
its  sublists. 
Term  -  A  list  containing 

(a)  A  constant 

(b)  A  sequence  of  pairs  corresponding  to  variables 
raised  to  powers. 

(c)  A  sequence  of  triples  -  each  corresponding  to  a 
function,  its  argument  and  a  powerj 

where 

(a')    The  constant  is  a  floating  point  number. 

(b')    The  first  element  of  each  of  the  pairs  is  a 

representation  of  a  variable,  e.g.,  a  Hollerith 
character.  The  second  element  is  a  floating  point 
number  representing  the  exponent  associated  with 
the  variable. 

(c')   The  first  element  of  each  of  the  triples  is  a  list 
name.   The  sublist  named  contains  the  expression 
which  is  the  argument  of  the  function. 
The  second  element  is  a  representation  of  the 
function  name. 

The  third  element  is  a  floating  point  number 
representing  the  exponent  associated  with  the 
function. 

In  Figure  1,  the  list  representation  of  the  expression  (l) 

is  diagrammed. 
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El 


TERM  1 


E2 


■TERM  ^ 


E 


TERM  1  --- 


TERM  2  --- 


1. 

2 

X             1 

2. 

0 

z           1 

1. 

0 

TERM  2 


1.5 


X 


3.1 


1.0 


E2 


SIN 


1.0 


E3 


SIN 


"W 


cW^ 


1.0 


I 


TERM  3  I 


1 

1. 

0         1 

X             1 

2, 

0           1 

E3 
I 

"term  n  \- 


TERM  4 


1 

1. 

0        1 

Y               I 

2. 

0 

e4 


TERM  5 


I 


i 


< — 


TERM  5  - 
TERM  b  - 


>* 

1. 

0        1 

X             1 

1, 

0 

TERM  6 


E5 


1.0 


Y 


1.0 


TERM  7  [---, 


E 


E5 


SIN 


1. 


TERM  7 


---^ 


1. 


Figure  1.   Representation  of  (1). 

Heavy  "black  mark  indicates  the  element  is  a 
sublist  pointer. 


-12- 


Ordering  within  Terms. 

Pairs  within  terms  appear  In  the  order  of  their 
associated  variables^  while  the  triples  which  appear 
within  terms  are  ordered  according  to  the  function  names. 
The  ordering  of  variables  is  on  the  numerical  value  of 
the  display  code  (on  the  66OO). 

A  Notational  Trick. 

Sometimes  a  user  may  want  to  work  with  an 
expression  such  as 

x.(x+y+z)200  . 

This  can  be  done  by  defining  a  function  f  with  the  property 
that  f(x)  =  X.   Then  the  above  expression  becomes 

Xf2°0(x+y+z) 

Which  is  of  the  type  easily  represented  by  triples.   In 
fact,  a  fimction  f  is  already  defined  within  SYMB66  and 
is  called  "1.*"  instead  of  "f".   The  appearance  of  "1.*" 
in  expressions  is  taken  into  account  in  the  relevant  parts 
of  the  code  such  as  the  output  and  derivative  routines. 
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B.      Input. 

The  most  convenient  method  of  introducing  an 
algebraic  expression  into  the  program  is  by  means  of 
input  data  cards.   This  method  is  discussed  here  first. 
An  alternative  method  is  then  discussed  for  creating 
expressions  during  the  execution  of  the  problem. 

The  expression  to  be  inputted  is  punched  onto 
cards  just  as  if  it  were  a  Fortran  expression  (with 
continuation  cards  if  necessary)  and  then  followed  by 
a  card  with  an  asterisk  in  column  1  (all  other  columns 
blank).   Every  expression  must  have  a  terminating  *  card. 
The  cards  are  arranged  as  data  in  the  order  the  programmer 
wishes  the  expressions  to  be  introduced  into  his  program. 
They  may  be  interspersed  with  other  data  cards. 

Then,  whenever  an  expression  is  to  be  inputted 
into  a  program,  the  following  statement  appears  in  the 
coding 

L  =   LISTPN(SL) 

where  both  L  and  SL  refer  to  the  same  expression,  as  if 
SL  appeared  in  an  EQUIVALENCE  statement  with  L. 

The  subroutine  LISTON  will  read  the  Fortran-like 
expression  on  the  next  data  cards  and  convert  it  into 
the  SLIP  list  discussed  in  section  A  so  that  the  expres- 
sion will  be  available  for  symbolic  computation.   For 
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every  expression  that  is  inputted,  a  separate  call  to 
LISTON  must  be  executed. 

To  create  an  expression  consisting  of  a  simple 
term  which  consists  of  a  coefficient  and  a  single  variable 
to  a  power,  i.e.,  of  the  form  C*  Z**PP  (where  C,  P  are 
floating  point  numbers)  and  call  the  resulting  expression 
Li;  the  sequence  of  statements  is 

LI  =  LIST(9) 

CALL  PUTLST(L1,C,1HZ,P,L1,L1) 
To  create  an  expression  (CZ*  FUNGT(Ll) *  * Q) 
consisting  of  a  simple  term  which  consists  of  a  coefficient 
C2  and  a  function  to  a  power  Q  with  an  argument  LI  (already 
our  expression)  and  call  it  L2,  the  statements  are 

L2  =  LIST(9) 

CALL  PUTLST(L2,C2,L1,5HFUNCT,Q,L2,L2) 
From  these  simple  expressions,  arbitrary 
expressions  usable  by  the  system  can  be  built  up  by  per- 
forming multiplications,  additions  and  substitutions. 
Warning:   An  expression  which  is  used  as  an  argument  of 
a  function  is  then  no  longer  available  to  the  user. 
If  it  is  necessary  to  keep  it,  use 

L8  =  LSSCPY(L7) 
l8  will  then  be  an  exact  duplicate  of  UJ . 
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C.      Output. 

The  bestj  simplest^  clearest  and  usual  way  to 
output  an  expression  is  to  use  the  routine  PRIPUT  which 
prints  expressions  as  card  images  of  Fortran  statements 
onto  tape  6.   A  typical  output  statement  is 

CALL  PRIPUT (6HRESULT, LA) 
where  the  first  argument  is  a  display  code  name  of  the 
expression  which  appears  to  the  left  of  the  equal  sign. 
The  second  argument  is  the  name  of  the  expression  to  be 
printed.   Obviously,  any  BCD  name  can  be  assigned  as  the 
first  argiunent. 

There  is  another  output  routine  which  may  be 
useful  in  iinusual  cases  such  as  removing  errors  or  display- 
ing lists  other  than  the  usual  algebraic  expressions.   The 
information  on  the  list  is  printed  in  octal  format  in  a 
vertical  array.   Such  a  format  is  hard  to  read  and  takes 
much  space.   For  these  reasons,  it  is  not  normally  used. 
The  calling  sequence  is 

CALL  PRLSTS(6HRESULT,LA) 
exactly  analogous  to  the  use  of  the  PRIPUT  routine.   It  is 
a  modified  form  of  the  SLIP  routine  of  the  same  name. 

Still  a  third  way  of  obtaining  output  is  to  use 
the  preprocessor  statement  OUTPUT  described  on  page  7 . 
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IV.     The  Algebraic  Computation  Package  of  Routines. 

A.     The  package  of  routines  is  conveniently  divided 
into  five  groups  according  to  the  types  of  operations 
they  perform. 

List  processing  routines:   These  are  the  SLIP  routines 
which  may  be  used  by  the  programmer  using  the  package. 

Expression  manipulation  routines:   Perform  algebraic 
and  arithmetic  operations  on  expressions  and  return 
expressions  as  final  results. 

Term  manipulation  routines:   Operate  on  terms  of  expressions, 
These  are  not  generally  needed  by  the  programmer. 

Operational  routines:   Perform  non-arithmetic  operations 
on  expressions  and  terms. 

All  other  routines   in  the  package  are  called  by  the  above 
routines.   A  list  of  these  other  routines  along  with  a 
brief  comment  about  each  is  given  near  the  end  of  this 
section.   The  programmer  need  not  concern  himself  with 
these  subalgebraic  routines. 
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1.     List-processing  routines  (SLIP) 

In  the  following,  list  operations  are  defined 

assuming  that  the  lists  are  algebraic  expressions. 

For  a  more  general  interpretation,  see  [2]. 

(i)     CALL  INITAS  (array  name,  dimension  of  the  array) 

This  call  must  appear  in  every  program  using  this 

package  of  subroutines.   The  subroutine  INITAS 

organizes  the  memory  allocated  by  the  user  for 

storage  of  expressions.   The  memory  requirement 

is  determined  by  the  user  in  a  DIMENSION  statement 

appearing  before  the  CALL. 

Example:  DIMENSI^^  ARR(IOOO) 

CALL  INITAS (ARR, 1000) 

(ii)   nom  =  LIST(9)   is  used  to  declare  the  name  of  an 

expression,  where  nom  is  the  name.   The  expression 
is  identically  zero  until  terms  are  added  to  it. 

(iii)      CALL  IRALST(nom) 

When  this  statement  is  executed  the  expression 
designated  by  nom  will  be  cleared  from  memory. 
All  further  references  to  the  expression  now  will 
be  meaningless  and  will  thus  cause  the  termination 
of  the  program. 

(iv)    CALL  MTLIST(nom) 

Makes  the  expression  nom  into  one  which  is 
identically  zero. 
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(v)    new  =  I^SGFY(nom) 

This  statement,  when  executed,  creates  an 
element -"by-element  copy  of  nom  which  is  then 
referred  to  by  new.   After  execution,  both 
new  and  nom  are  the  same  algebraic  expressions 
although  they  occupy  different  locations  in 
memory. 

(vi)    ival  =  LSTEQL(jga,ib) 

LSTEQL  compares  the  expression  ^  with  the 
expression  ib_     element-by-element.   If  the  two 
expressions  are  identical,  but  perhaps  in  different 
parts  of  computer  memory,  then  the  value  of  LSTEQL 
is  zero.   Otherwise  it  is  non-zero. 

(vii)   ival  =  LISTiyiT(nom) 

LISTMT  tests  to  see  if  the  expression  named  nom 
is  currently  zero.   If  the  expression  is  zero, 
ival  =  0.   Otherwise  ival  y  0. 
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2.     Expression  Manipulation  Routines. 

The  variables  used  in  the  given  descriptions  are 
designated  by  the  following  symbols  to  indicate  the  type 
of  variable  required  or  produced  by  the  routine.   The 
prefix  "x"  merely  indicates  the  floating  point  named 
equivalent  of  the  variable  without  the  prefix  (no 
conversion).   All  the  routines  may  be  used  as  functions 
or  subroutines. 

(1)  char    Hollerith  variable 

(2)  fpo     floating  point  number 

(3)  lad        the  calling  argument,  which  is  the  name 

of  an  expression  to  be  operated  upon,  is 
not  preserved,  i.e.  is  destroyed. 

(4)  iap        the  calling  argument,  which  is  the 

name  of  the  expression  to  be  operated  upon, 
is  preserved  throughout  operation. 

(5)  ibp 

(6)  ico     calling  expression  which  is  operated  upon 

by  the  routine  to  become  resulting  expression 

(7)  icr     resultant  expression^   created  by  routine 

if  not  already  designated  by  programmer. 
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(1)  xicQ  =  ADD(  iad,  ico)  adds  iad  to  ico 

(2)  xicr  =  DVSUM( ^ap, char, ^cr )  divides  jap  by  char 

and  adds  result  to  icr 

(3)  x^co  ^   EVALUE(  ^co^v-,  jGon-.  ,  .  . . ,  v^^con^)  evaluates 

ico  for  the  variable  v.  equal  to  the 
constant  con.  (1  j<  i  j<  6). 
The  following  functions  are  automatically 
evaluated,  when  their  arguments  are  constant^ 
sin,  cos,  exp,  log,  tan,  arctan  and  tanh. 

(4)  xicr  =   GETCg(E(char,fpo,  jap,  ^cr)  places  the 

coefficient  of  char  **  fpo  on  jcr. 

(5)  xicr  =   H1TENT( ^ap,char, jbp, jcr)  replaces  jap  for 

char  wherever  it  appears  in  ibp  and  adds 
results  to  icr.   icr  results  in  factored 
form  and  no  replacement  is  done  in  function 
arguments. 

(6)  icr  =   INSBST(  jap,  char,  ^bp,icr)  substitues  jap  for 

char  wherever  it  appears  in  £bp   and  adds 
result  to  icr.  icr   is  not  in  factored  form. 

(7)  -gcr  =   I]NfSUBT(iap,char,  ibp,  icr,Vj^,Pj^,Vg,Pg, .. . 

(up  to  V  ,p  ))  substitutes  as  in  (6)  but 
eliminates  all  terms  containing  the  variable 
V.  **    d,  where  v.  is  a  Hollerith  variable  and 
d  is  a  floating  point  exponent  ^  p.. 
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(8)  icr  =    INTENT( iap, char, £bp,  icr)  replaces  as  in  (5) 

but  includes  function  arguments. 

(9)  mc  =   LDIF( iapjChar)  differentiates  jap  with 

respect  to  char  and  creates  the  resulting 

expression  mc_.   The  derivatives  for  the 

following  fiinctions  have  been  built  into 

the  package.   The  user  may  add  other  functions 

and  their  derivatives  to  the  system  for  his 

own  purposes  by  using  the  routine  LOOKUP 

(see  page  25). 

sin     cos     exp     log     tan 

cot     sec     CSC     arcsin  arccos 

arctan  artanh  arsech  arccot  arcsec 

arccsc   sinh    cosh    tanh    coth 

sech    csch    arsinh  arcosh  arcoth  arcsch 

(10)  £co   -   NUMFY(ico,fpo)  multiplies  ico  by  fpo. 

(11)  xicr  =  SBST( jap, char, £bp, icr)  substitutes  as  in  (6) 

but  does  not  consider  f\inction  arguments. 

(12)  xicr  =  Sg'LVE(  ^apjChar,  icr)  assumes  ^ap  to  equal  zero 

(not  identically),  and  that  char  appears 
linearly  in  exactly  one  term  of  iap.   S{2lLVE 
solves  for  char  and  places  result  on  Zcr. 

(13)  xico  =  SUB( ^ad, ico)  subtracts  ^ad  from  ico. 

(14)  xicr  =   SUBT(  jap, char,  ibp,  ^cr, v-^^,p.j^,  .  .  .  (up  to  v  ,p  )) 

substitutes  as  in  (7)  but  does  not  consider 
function  arguments. 

-22- 


(15)  xicr  =   SUMPY(  jap,  ibp,  £ct)   multiples  ^ap  by  ibp 

and  adds  results  to  icr. 

(16)  x^co  =   TRCAL( ico, tchar,Ghar)  finds  smallest  power 

pwr  of  char  in  the  terms  in  which  the 
variable  tchar  is  contained,  and  then 
truncates  all  terms  with  tchar  to  powers 
greater  than  pwr  in  ico. 

(17)  x.gco  =   TRIINC(  ^co,charjfpo)  eliminates  all  terms  in 

ico  whose  power  of  char  is  greater  than  fpo, 
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5.      Term  Manipulation  Routines. 

(1)  vai,    =       CPYTRM(  £a,  £b)      If  £a_   is  not  empty,  one  term 

is  removed  from  expression  £a_   to  expression 
it)  and  vai  =  1.   If  jga  is  empty,  va£  =/  1. 

(2)  x£      =       DVTRM(  i,char)   divides  term  _i  by  char  where 

£_   is  the  original  as  well  as  the  resulting 

term  x£. 
(5)    xl   "   SRTRM(^)  arranges  the  variables  of  the 

term  _Z   (equiv.  x£)   in  lexicographical  order. 

Then  follows  the  variables  with  the  fvmctions 

which  are  also  ordered. 
(4)    xl   ^   TRMMPY(m,n,^)  multiplies  term  m  by  term  n 

to  make  term  £,   which  is  equivalent  to  xi. 

If  _£   has  not  been  made  a  list  name,  TRMMPY 

also  makes  £   a  list  name. 
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4.      Operational  Routines. 

(1)  xia  =   BREAK(  £,  char,  ib,  ia  )  separates  the  given 

expression  _£   into  expressions  £a_   and  Jb. 
Those  terms  which  contain  char  are  placed 
on  £b_,    and  those  terms  which  do  not,  are 
placed  on  £a_   (equivalent  to  xia ) .   If  _£a_ 
and  £b_   have  not  previously  been  designated 
as  expressions,  BREAK  declares  them. 

(2)  £a     =        LEXIC(^(xia)  collects  the  terms  of  the 

expression  _£a_  (equiv.  to  xia )  and  also 
places  them  in  a  conventional  lexicographical 
order. 

(3)  ia     =        LIST{2'lN[(xia)  creates  a  SYMB66  expression  from 

a  data  card  (see  page  l4). 

(4)  CALL  L^pKUP(-l. ,fname, Mer)  allows  the  user  to  make 

available  to  LDIF,  the  differentiation  routine, 
derivatives  of  functions  used  in  his  program, 
which  are  not  already  built  into  the  package. 
(For  a  list  of  fxinctions  whose  derivatives 
have  been  built  in,  see  page  22).   fname  is 
the  Hollerith  name  of  the  function  for  which 
the  derivative  is  being  specified.   ider  is 
the  name  of  the  "expression"  created  for 
the  derivative  of  fname.   Zero  is  always 
given  in  place  of  the  argument  of  fname  in 
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"expression"  ^der .   For  example  to  specify 
the  derivative  of  tan  (arg)I 

1  LL  =  LIST(9) 

2  CALL  PUTLST(0,LL,1-^0,3HSEG,2.,LL,LL) 

3  CALL  L?f0<UP(-l.,3HTAN,LL) 

Note  in  statement  2,  zero  is  the  first 
argument  to  PUTLST.   This  is  required 
whenever  PUTLST  is  used  to  create 
derivative  expressions. 

(5)  fpo  =   PgiWER(ia,char)  obtains  the  power  fpo  of 

the  variable  or  function  char  where  it 
first  appears  in  the  expression  ia. 

(6)  CALL  PRLSTS (Hollerith  name,  Ja)  prints  out  any  list 

designated  _ia  in  octal  format  (see  page  l6). 

(7)  xia  =   PRIPUT  (Hollerith  name,  ia_)   prints  out 

any  valid  expression  in  Fortran  notation 
(see  page  l6) . 

(8)  xia  =   PUTLST(ia,a,b,c, . . . , ia, ^a)  adds  terms  to 

the  expression  i^.      The  last  term  given  in 
the  calling  sequence  is  followed  by  a  double 
ia,  and  each  term  is  followed  by  lo^.      See 
Section  III  on  Notation. 
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(9)  xia  =   SMPL( £a)      orders  the  terms  of  the  expression 

ia  within  themselves  and  then  collects  the 
terms  and  places  them  in  lexiographical 
order. 

(10)  vai  =    TSTG0J(ia,Gon)   tests  if  expression  £a_   is 

equal  to  a  constant.   If  _ia_  is  not  equal 
to  a  constant  vai  =0.   If  ia  is  equal  to  a 
constant,  including  the  zero  case,  v&£  =  1., 
and  con  is  the  constant. 
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5.     Sub algebraic  Routines. 

(a)    These  routines  are  SLIP  routines  and  descriptions 
for  them  can  be  found  in  [2]. 

LNKR  B^  NAMTST  STRIND  SEQRDR 

P0PT(2fP  AND  LNKL  P(2fPB^  DELETE 

ID  SEQLR  SEQLL  LCNTR  MAPgfV 

SETDIR  NUCELL  C(^T  RCELL  NEWT0P 

STRDIR  INHALT  NEWB^^I^  L(2fCT  SETIND 

T0t  EQUAL 

Some  of  these  routines  were  carefully  rewritten  in 
ASCENT  (GDC  6600  assembly  language)  which  saved  a  little 
machine  time.   A  more  powerful  technique  has  proven  to  be 
the  overwriting  of  transfers  to  short  closed  routines  with 
an  equivalent  open  code.   Doing  this  overwriting  at  execu- 
tion time  has  proven  to  be  both  simple  and  efficient, 
yielding  significant  savings  in  central  processor  time. 

(b)    The  first  three  below  are  list  processing  type 
routines : 

(1)  LSCMPR  compares  the  elements  of  two  modified 
expressions  for  equality. 

(2)  P$2('PMID  deletes  the  list  element  which  is  currently 
being  pointed  to  by  the  sequence  reader  and  the 
value  assigned  to  the  function  is  the  value  of 

the  element  deleted. 
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(3)    MANY(nom,a,b,Cj . . . )   up  to  59  variables  are  appended 
to  the  bottom  of  the  list  nom  in  the  order  the 
variables  are  given,   nom  must  already  be  the  name 
of  a  list. 

The  routines  below  are  used  by  other  routines  in  the  package. 

{H)  HREVAL  is  used  by  EVALUE 

(5)  RAISE   is  also  used  by  EVALUE;  RAISE  performs 
exponentiation 

(6)  F<ilfFX         evaluates  a  function  with  a  numerical 
argument. 

(7)  EXTRAC  determines  if  a  given  variable  is  present 
in  a  specified  term  ,     and  if  so,  to  what  power. 

(8)  KEBLA,  FILL,  NIPF   are  used  by  PRIPUT 

(9)  KgfKE  is  used  by  PUTLST 

(10)  INLIST  is  an  input  routine  used  by  the  preprocessor 
QUEST 

(11)  LCPYDR  is  used  by  LDIF  to  form  the  deriative  of 
a  fimction 

(12)  SAME,  m0N0       are  used  to  make  Fortran  fixed-floating 
point  conversions 

(15)    ADD IN   is  used  in  lexicographical  ordering 
(l4)    ERRFUT  is  used  to  print  out  error  messages  and 
abort  execution. 
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B.     An  Illustration. 

Suppose  sin(x)  is  to  be  evaluated  in  its  power  series 
expansion  for  ten  terms  and  x  =  2. 

sin(x)  =  x-yr  +  -^--.-+  -^^   . 

The  coding  for  a  complete  program  would  be  I 
PR{2^GRAM  SINGH  (INPUT,  OUTPUT) 

cgfimgfN  avsl 

DIMENSIi^t^  ARR(IOOO) 

C  {2IMM0^ /SYMS  YS /S  PQR  ( 1 0  00 ) 

CALL   INITAS(ARR,1000) 

NSIG  =    -1. 

NFAC   =   1. 

K  =   19 

LSIN   =   LIST(9) 

D^  1      J=1,K,2 

NSIG  =   NSIG  *     (-1) 

IF(J.NE.1)NFAC  =    (NFAC+2)*(NFAG+1)   *    NFAC 

FJ  =   J 

FF  =   NSIG*     (1./PLPAT(NFAC) 

LSN  =   LIST (9) 

GALL  PUTLST(LSN,FF,1HX,FJ,LSN,LSN) 
1    CALL  ADD (LSN J LSIN) 

CALL   PRIPUT ( 6 HS IN ( X ) , LS IN ) 
CALL  EVALUE(LSIN,1HX,2.) 
CALL   PRIPUT(6HSIN(2.  ),I^IN) 
CALL  EXIT 
END 
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V.      Preprocessor  (QUEST) 

A  simple  problem  oriented  language  and  preprocessor 
called  QUEST  (QUick  and  Easy  Symb66  Translator)  has  been 
designed  to  facilitate  the  use  of  Symb66.   This  preprocessor 
translates  into  Fortran.   Since  the  target  language  is 
Fortran,  QUEST  statements  can  be  freely  intermingled  with 
Fortran  statements.   By  coding  in  QUEST  the  programmer  is 
relieved  of  the  task  of  learning  the  internal  structure  of 
Symb66  and  the  various  Fortran  subroutine  calls.   For 
symbolic  differentiation  or  substitution,  the  use  of  QUEST 
may  be  especially  convenient. 

The  rules  for  the  form  of  QUEST  statements  are 
given  below  followed  by  information  on  error  messages  and 
instructions  for  programming  in  QUEST.   (A  list  of  QUEST 
statements  with  their  corresponding  SYMB66  subroutine  names 
and  a  sample  QUEST  program  is  given  in  Section  II  of  this 
report  (page  5 ) . ) 

A.     Format  of  statements. 

(1)  Statement  ninnbers  and  continuation  cards  are  written 
in  the  usual  Fortran  manner,  with  up  to  19  continuation 
cards  allowed.   Statements  begin  in  column  7. 

(2)  Blanks  are  ignored. 

(3)  No  variable  identifier  longer  than  ten  characters  is 
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accepted. 

(4)  All  fixed  point  numbers  are  automatically  converted 
to  floating  point  numbers  in  QUEST  statements. 

(5)  QUEST  statements  may  be  interspersed  with  Fortran 
statements,  including  Symb66  subroutine  calls. 

(6)  Every  variable  identifier  (char,  con,  name,  pur, 
var,  term,  etc.)  must  be  followed  by  a  comma  (,)  or  an 
equal  sign  (=),  except  where  the  variable  identifier 
teraiinates  an  executable  statement.   (There  is  no  terminal 
piinctuation. ) 

(7)  An  X  in  column  1  indicates  the  Fortran-like 
arithmetic  expression  is  a  symbolic  expression. 

(8)  An  S  in  column  1  indicates  that  the  desired 
arithmetic  operations  are  symbolic. 

B.     Error  messages. 

Unacceptable  QUEST  statements  are  flagged  with  an 
error  message.   The  errors  recognized  are! 

** comma**         comma  missing  after  an  identifier 

** final  punctuation**   statement  not  correctly  terminated 

** constant**      constant  incorrect  or  not  found 

when  expected 
**instinjiction**   use  of  QUEST  vocabulary  incorrect 
** equal  sign**    equal  sign  missing  when  expected. 
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A  listing  of  the  QUEST  program  is  always  given  to 
the  programmer  when  errors  are  found;,  and  the  number  of 
errors  appears  at  the  end  of  the  listing.   Also  the  message 
"QUEST  DR^S  J^"  appears  in  the  user's  dayfile,  and  the 
user's  entire  job  is  aborted,  so  that  no  further  processing 
takes  place. 

C.      Instructions  and  Information  for  Programming  in  QUEST 

(1)  TAPE  2  must  not  be  used  by  the  programmer,  and 
the  variable  names  ABYSS,  AVSL,  SPQR,  and  SYMSYS. 
(QUEST  uses  the  file  name  TAPE  2  and  the  variable  names.) 

(2)  A  QUEST  program  begins  and  ends  like  a  Fortran 
program.   However  the  user  must  include  the  "RESERVE" 
statement  followed  by  the  total  memory  locations  he  has 
estimated  for  manipulating  the  algebraic  expressions  in 

his  program.   The  RESERVE  statement  follows  the  header  card. 

(3)  When  the  QUEST  translator  has  processed  a  program 
successfully,  the  message  "QUEST  C^PLETES  J{2te"  appears 
in  the  user's  dayfile.   Subsequent  control  cards  are  then 
executed. 

(4)  All  expressions  appearing  in  the  executable  group 
of  statements  must  have  previously  been  declared  (using 
EXP,  X  notation)  or  have  been  created  in  execution  of  a 
previous  instruction.   This  latter  case  includes  the 
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expression  identifier  following  AND  PLAGE  0   (AP^ . 
If  these  identifiers  have  not  been  previously  declared 
to  be  expressions,  the  identifiers  will  be  made  into  the 
desired  expressions  during  execution  of  the  statement, 
and  will  then  be  available  for  use. 

(5)  All  QUEST  expressions  are  preserved  throughout 
execution  of  program  except  when  specifically  changed 
or  destroyed.  Exceptions  to  this  statement  are  noted 
in  Section  II. 

(6)  QUEST  abbreviations  may  be  freely  substituted  for 
the  longer  forms  of  the  statement.   See  Section  II, 

D.     QUEST  statements. 

Preliminary  comments  on  notation. 

(1)  All  abbreviations  are  either  underlined  or  appear 
in  parentheses  following  the  appropriate  phrase.   As  a 
rule  the  longer  commands  are  abbreviated  to  the  first  three 
letters  of  the  word,  and  phrases  of  three  words  are  abbrevi- 
ated to  a  mnemonic  consisting  of  the  first  letter  of  each 
of  the  three  words. 

(2)  Several  statements  in  both  the  executable  group 
and  the  declarative  group  may  be  written  as  a  standard 
type  Fortran  statement  using  the  customary  arithmetic 
operators.   QUEST  recognizes  that  these  statements  are 
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symbolic  (rather  than  the  usual  Fortran)  by  an  S  (indicat- 
ing operation)  or  X  (indicating  declaration)  placed  in 
column  1  by  the  programmer. 

(3)  In  the  examples  below  '^a'^,  '  ib '  ;,  and  '  ^c '  are  all 
expression  names  and  are  variable.   The  SYMB66  or  SLIP 
equivalent  for  each  QUEST  statement  follows  the  semicolon.* 

(4)  The  boldface  characters  are  the  elements  of  the 
QUEST  language. 

E.      QUEST  instructions. 

Declarative-type  statements 

(1)  To  create  a  symbolic  expression: 
EXPRESS  I  glH  ia  =  Fortran  type  expression 
(use  of  "/'  operator  prohibited) 

or  X  in  coliomn  1  and  merely  the  expression  ia 
beginning  in  column  7j 
LIST0^ 

(2)  To  designate  the  name  of  an  expression,  whose 
value  is  imderstood  to  be  zero  until  other  terms  are  added.' 

NAME  Za,m,ic,...    (up  to  5  var); 
LIST 


The  semi-colon  is  not  part  of  QUEST. 
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(3)  To  add  a  derivative  for  a  function  not  built  into 
the  SYMB66  system 

THE  DERIVATIVE  FgfR   (TDF)  function  name,  is 
£a   =  expression  written  with  zero  arguments 
(no  "/"s); 

(4)  To  allocate  memory  locations  for  symbolic 
manipulation! 

RESERVE  integer;       C{2{^^  AVSL  and 

cpMM(2(l^/SyMSYS/SPQR(1000)  and 

DIMENSIpi^  ABYSS  (n)  and  CALL  INITAS  (ABYSS, n) 

F.     Executable  type  statements. 
a.     Arithmetic  operations 

(1)  To  add: 

ADD  £a,    T^  £b',    ADD  -  expression  £a   is  destroyed. 

(2)  To  subtract: 

SUBTRACT  ia,FR^  ib;  SUB  -  £a   is  destroyed. 

(3)  To  divide  an  expression  by  a  character: 

DIVIDE  ia,  BY  char,  AND  PLACE  0    {AFgi)    ib;  DVSUM. 

(4)  To  multiply  an  expression  by  another  expression: 
MULTIPLY  ia,  BY  ib,  AP^  £c;    SUMPY. 

(5)  To  multiply  an  expression  by  a  constant 

MUL  £a,    BY  con;  NUMPY  -  £a   is  the  original  and 
resulting  expression. 
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The  above  instructions  (1-5)  niay  all  be  written 

in  the  following  manner  with  an  S  in  column  1^ 

ib  =  ib  +  ia  ADD^  SUB 

m   ^  Ha  /   char  DIV 

!tc    =  Ha  *     ^b  MUL 

ia  =  ia  *  con  MUL 

(6)  To  differentiate  an  expression  with  respect  to 

a  char;  check  that  derivative  is  in  system  (see  page  22) 
or  use  the  TDF  statement! 

DIFFERENTIATE  £a ,    WITH  RESPECT  T^  (WRT)  char, 

AP(2f  m;    LDIF. 

(7)  To  evaluate  a  symbolic  expression  by  giving  a  value 
to  one  or  more  variables! 

EVALUATE  i,a ,    F^R  char  =  no,  char  =  no,... 
(up  to  5  char  =  no);  EVALUE  -  ia    is  the 
original  expression  as  well  as  the  evaluated 
expression   (for  functions  automatically 
evaluated  see  page  21). 

(8)  To  obtain  the  coefficient  of  a  character! 

GET  C^FFICIENT  g^    (GC^)  char  **  no,  IN  Ha,    AP^  m', 
GETCi2te  (if  **    no  is  omitted  the  exponent  is 
assumed  to  be  1 ) . 

(9)  To  obtain  the  exponent  of  a  character! 
GET  P(2IWER  pF  (GPi2f)  char,  IN  Ha,    AFgi   var; 
PpWER   (var  is  a  variable). 
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(10)  To  eliminate  terms  in  an  expression  that  contain 
char  to  a  power  greater  than  no. 

TRUNCATE  ia,  F0R  char  **no;  TRUNC  ia  is  the 
original  and  resulting  expression;  see  (8) 
for  explanation  of  char  **no. 

(11)  To  solve  for  char  assioming  the  given  expression  is 
equal  to  zero,  and  that  char  appears  linearly  in  only  one 
term  of  the  given  expression: 

S^VE  ia,  FgfR   char,  AFgf   ^t;  S0LVE  . 

The  next  six  statements  (12-17)  all  pertain  to  substitution. 
To  clarify  the  differences  among  the  different  statements, 
bear  in  mind  that  the  word  REPLACE,  as  it  is  used  here, 
means  to  merely  insert  an  expression  into  the  original  in 
place  of  a  designated  character,  but  not  to  carry  out  any 
multiplication.   While  on  the  other  hand,  SUBSTITUTE  means 
to  insert  an  expression  for  a  specified  character  and  to 
carry  out  all  necessary  multiplications  in  order  to  produce 
a  simplified  version  of  the  original  expression. 

(12)  To  substitute  an  expression  £a_   for  char  appearing 
in  ib,  where  appearances  of  char  in  function  arguments 
is  ignored: 

SUBSTITUTE  £a,    FgfR   char,  IN  ib,  AFg!   ic,  EXCLUDING 
FUNCTI^  ARGUMENTS  (EFA);  SBST  . 

(13)  To  substitute,  as  in  (12),  but  also  make 
substitutions  for  char  appearing  in  function  arguments! 
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SUB  ia,  Fi2(R  char,  IN  £b,    A.P{2^  ic;  INSBST. 

(14)  To  substitute  as  in  (12),  but  to  eliminate  all 
terms  containing  var  to  a  power  greater  than  the  specified 
pwr: 

SUB  £a,  Fpi?  char,  IN  £b,  APgf  Ic ,  EFA  AND  TRUNCATE 
0T  (AT^  var,  pwr,  var,  pwr,  ...  (up  to  5  pairs); 
SUBT  . 

(15)  To  substitute  as  in  (l4)  but  to  include  function 
arguments  as  in  (15)* 

SUB  ia,  Fpfe  char,  IN  £b,  kV^  ic-  AT{/ var,pwr, . . . 
(up  to  5  pairs);  INSUBT  . 

(16)  To  replace  char  whenever  it  appears  in  _|al 
REPLACE  char,  IN  ia,  WITH  ib,  kV(^   ic;  INTENT  . 

(17)  To  replace  char  in  ia,  but  excluding  appearances 
of  char  in  f\inction  arguments  T 

REP  char,  IN  ia,  WITH  ib,  AP0'  ic,  EFA;  HITENT  . 

b.     Non-arithmetic  expressions. 

(1)  To  print  out  one  to  five  expressions! 
0tJTPUT  ia,  ib,  ...  ;  PRIPUT  . 

(2)  To  eliminate  all  terms  of  an  expression,  so  that 
the  expression  is  now  equal  to  zero  I 

EMPTY  ia, ib,...  (up  to  5);  MTLIST  . 

(3)  To  erase  from  memory  an  entire  expression,  so  that 
even  the  name  can  no  longer  be  referenced! 
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DESTRi2^  Za,    £h ,    ...  (up  to  5);  IRALST  . 

(4)  To  make  a  copy  Zb_   of  the  expression  £a,    where  £a   =  ib 
C(2fFY  ia,  ig(   ib;  LSSCPY  . 

(5)  To  simplify  an  expression,  including  its  fiinction 
arguments,  by  collecting  terms: 

SIMPLIFY  £a,    ih ,    ...     (up  to  5);  SMPL  . 

(6)  To  order  the  terms  of  an  expression! 
S^T  £a,    ib,  ...  (up  to  5);  LEXlCgl  . 

This  instruction  is  rarely  used  since  lexicographical  order 
is  preserved  by  the  SYMB66  system,  and  the  user  of  QUEST 
ordinarily  need  not  concern  himself  with  this  aspect. 

Tests. 

(1)  Constant  test  -  if  the  expression  £a   is  equal  to  a 
constant,  the  value  of  the  constant  is  stored  in  the  Fortran 
location  varj  in  either  case  control  may  be  transferred  to 
the  next  statement j  or  if  a  statement  no.  appears,  control 

is  transferred  to  that  statement  if  £a   is  equal  to  a  constant^ 

if  ELSE  follows  the  statement  no.  control  is  transferred 

to  the  number  following  ELSE  if  £a   is  not  equal  to  a  constant! 

IF  £a,    C{^STANT  var  =  cgfN 

(optional)  statement  no. 

(opt.),  ELSE  statement  no.;  TSTG0^  . 

(2)  Zero  test  -  if  £a   is  equal  to  zero,  control  is 
transferred  to  statement  no.  following  7,ER0',    if  £a    is 


-40- 


not  equal  to  zero,  control  is  either  transferred  to  next 
statement  or  statement  following  ELSE  if  ELSE  appears." 

IF  Ha,    ZER^  statement  number  (optional),  ELSE 

statement  no.;  LISTMT  . 
(3)     Equality  between  two  expressions  -  if  ^a  is  an 
exact  copy  of  ib  control  is  transferred  to  given  statement 
number;  if  ia  is  unequal  to  ib,  control  goes  to  next 
statement,  or  statement  following  optional  ELSE; 

IF  Za,    EQUAL  m ,    statement  no.  (optional), 

ELSE  statement  no.;  LSTEQL  . 

G.      Examples  of  QUEST  Statements  . 

The  following  merely  demonstrates  valid  QUEST 
statements.   Although  the  program  was  compiled,  it  is 
not  intended  to  be  mathematically  meaningful. 
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Rfg^Sv?  MAJg'^WPUT^ 'OUTPUT) 

EXPLA»2+7.5»X  +A*Y»«2«Z 
EXPRESSION  LB  =    15»X»»2 
COPYLA,  TO  LC 
ADD  LA*  TO  LB 
CALL  PRIPUT(2hLB.LB) 
DIF  LCtWRTX.APOLD 
CALL  PRIPUT(2HLD»LD) 
DIVIDE  LCtBY  Y»APOLE 
CALL  PRIPUT(2HLE»LE) 

DES  LE 

EVALUATELD^FOR  X  -  7.5 

CALL  PRIPUT  (2HLD.LD) 

EMPTY  LB 

DESTROY  LCtLB.LO 

EXP  MA«3»X»Y*»2»Z»«3+W-A»B»«2. 

OUTPUTMA 

COPYMA.TO  MB 

1  EVA  MA,FOR  X=2.»B«10.   W-16 
4  OUTPUT  MA»MB 

2  DIFFERENTIATE  MB»WITH  RESPECT  TO  1*    AND  PLACE  ON  MC 
OUT  MC 

3  NAME  MD»ME 
OUT  MD  ,  ME 

EXP  MF.4»SlN(pn»C0S(X4-2»Pn  +3»PI 
OUT  MF 
6  REPLACE  W*tN  MB»WITH  MFtAPO  MDtEXCLUDING  FUNCTION  ARGUMENTS 
OUT  MO 

REP  PI .IN  MDt  WITH  MC .  ANO  PLACE  ON  ME 
OUT  ME 
EXP  Ll«D*»3 
OUT  LI 

EXP  LA«A+B-^C 
ADD  L4.  To  L2 
OUT  L2 

SUBL2»  FOR  D»  IN  LI.  APO  L3»  EFA 
OUT  L3 

SUB  L2»  FOR  D»  IN  LI •  APOL  A»  ATO  A»3«»6»3 
OUT  L4 


+  20*A 


X 

LA  «  X»»6 

X 

LB  =  15»A»X»»2 

s 

LB  =  LA  +  LB 

OUT  LB 

X 

LA  =  12»A»X 

s 

LB  =  LB  -  LA 

OUT  LB 

COPY  LB.  TO  LD 

OUT  LD 

s 

LA  s  LB/A 

OUT  LA 

s 

LA  =  LA  ♦  2 

OUT  LA 

s 

LB  »LA*.5 

OUT  LB 

s 

LC  =  LB»LA 
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OUT  LC 

IF  LD»EOUAL  LC.  10 

OUT  LD 

OUT  LC 

10  DESTROY  LDtLA.LB 

GCO  X»»2»IN  LC.  APO  LA 

OUT  LA 

GET  POWER  OF  A. IN  LA.  APO  PWR 

PRINT  11. PWR 

11  F0RMAT(1H1.«   PWR  =  ♦F10.5//) 
LPOW  =  1, 

CALL  NUMPY(LPOW.PWR) 
OUT  LPOW 
DES  LA 

DIFFERENTIATE  LC.WRT  X.  APO  LA 
OUT  LA 

LB  »  C0S(2»B  +  Y) 
OUTPUT  LB 

MULTIPLY  LA.  BY  LB. APO  LD 
OUT  LD 

EVALUATE  LD.  FOR  B  «  1. 
OUT  LD 

EVA   LD.FOR  Y=-2. 
OUT  LD 

COPY  LD.  TO  LE 
EVA  LD.  FOR  X»  1..A  «  2. 
OUT  LD 
EWPTY  LD 
ADO  LE,  TO  LD 
OUT  LD 

REPLACE  A  »  IN  LC  WITH  LB.  APO  LE 
OUT  LE 
LF  a  Y 

REP  B,  IN  LE.  WITH  LF.  APO  LG.EFA 
OUT  L6 

DES  LE.  LG*  LF 
K  »  10 

NAME  LEXtMEX 
NS16  «  -1 
NFAC  »  1 
D0710   J«1.IC.2 
NST6  «  NSIG»(-1) 

IF  (J  .NE.  1)   NFAC  »  (NFAC+2)*(NFAC+1)»NFAC 
FJ  *  J 
FAC  a  NFAC 
FF  -  NSIG»(1«/FAC) 

CALL  PUTLSTILEX.FF.IHY.FJ.LEX.LEX) 
710  OUT  LEX 

NFAC  «  1 
DO  20   Jal.K 
FJ  »  J  -  1 

IF  (FJ  .NE.  0)   GO  TO  19 
CALL  PUTLST(MEX.1..MEX»MEX) 
GO  TO  20 
19  CONTINUE 

IF  (J  ,NE.  1)   NFAC  «  NFAC»J 
FF  «  l./NFAC 

^5 


CALL  PUTLST(MEX.PF,lHZ»FJ»«EX.MEy) 
20  OUT  MEX 

SUB  LEX, FOR  Z,  IN  MEX.APO  LEXSIN»EFA 
OUT  LEXSIN 

END 
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person  acting  on  behalf  of  the  CoianiBslon: 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this 
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